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We extend the quantitative string evolution model of Martins and 
Shellard to superconducting strings by introducing a simple toy model 
for the evolution of the currents. This is based on the dynamics of a 
‘superconducting correlation length’. We derive the relevant evolution 
equations and discuss the importance of plasma effects. Some conse¬ 
quences of our results are also suggested. 
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I. INTRODUCTION 

It is well known that cosmic strings can in some circumstances (typically when 
the electromagnetic gauge invariance is broken inside the string) behave as ‘super¬ 
conducting wires’ carrying large currents and charges—up to the order of the string 
mass scale in appropriate units |jj. The charge carriers can be either bosons or 
fermions (see [^| for a review). The former type occurs when it becomes energet¬ 
ically favourable for a charged Higgs field to have a non-zero vacuum expectation 
value in the string core; the latter happens when there are fermion zero modes 
moving in the string field. 

If superconducting strings carry currents, they must also carry charges of similar 
magnitude. This includes not only charges trapped at formation by the Kibble 
mechanism but also the ones due to string inter-commuting between regions of the 
string network with different currents. However, arbitrarily large currents or charges 
are not allowed—there are critical values beyond which the current saturates and 
charge carriers can leave the string. 

Even though the overwhelming majority of the work done on cosmic strings so 
far was concerned with the structureless Goto-Nambu strings, there are a number 
of cosmological scenarios in which superconducting strings play a crucial part— 
notably the so-called ‘vorton problem’ ||— so that their study is highly relevant. 

In this paper we take a first step in this direction by generalizing the quantitative 
string evolution model of Martins and Shellard 0 |^] . After a brief review of this 
model, we extend it by developing a simple toy model for the evolution of the 
currents on the long strings, based on the dynamics of a ‘superconducting correlation 
length’. The different contributing mechanisms are discussed, and in particular we 
consider the importance of plasma effects. Finally, we briefly mention some of 
the implications of our model—a more detailed discussion is left for a forthcoming 
publication! |J. 
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C. J.A.P.Martins @ damtp.cam.ac.uk 
^Electronic address: E.P.S.Shellard@damtp.cam.ac.uk 
Submitted to Phys. Lett. B. 
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Throughout this paper we will use fundamental units in which Ti = c = ks = 

Gm 2 pi = 1 . 


II. EVOLUTION OF THE CURRENTS 


Analytic string evolution methods describe the string network by a small number 
of macroscopic quantities whose evolution equations are derived from the micro¬ 
scopic string equations of motion. The first such model providing a quantitative 
picture of the complete evolution of a string network (and the corresponding loop 
population) has been recently developed by Martins and Shellard [j|j|] , and has two 
such quantities, the long-string correlation length = p/L 2 (p being the string 
mass per unit length) and the string RMS velocity, v 2 = (x 2 ). It also includes two 
‘phenomenological’ terms, a ‘loop chopping efficiency’ parameter 0 < c < 1/2 and 
a ‘small-scale structure parameter’ 0 < k < 1. Their evolution equations are 

2—= 2HL(1 + v^) + v^—+ cvoo , (2.1) 


dv 

dt 


(1-u 2 ) 



2H+ £ f 


( 2 . 2 ) 


these are sufficient to quantitatively describe the large-scale properties of a cosmic 
string network. 

Note that the explicit form of the ‘friction lengthscale’ if will depend on the type 
of current involved. For the case of a neutral current, one expects Aharonov-Bohm 
scattering ]l0| to be the dominant effect, and consequently we have [|| 

e ' = W>- i 2 - 3 ' 


where T is the background temperature and (3 is a numerical factor related to the 
number of particle species interacting with the string. We assume that this is the 
case for the time being, but the more general case where plasma effects can be 
important will be discussed towards the end of this paper. 

We now extend it setting up a ‘toy model’ for the evolution of the currents. 
One assumes that there is a ‘superconducting correlation length’, denoted £, which 
measures the scale over which one has coherent current and charge densities on the 
strings. Associated with this we can define N to be the number of uncorrelated 
current regions (in the long-string network) in a co-moving volume V, and it is then 
fairly straightforward to see how the dynamics of the string network affects N and 
obtain an evolution equation for it. 

Firstly, we expect that in a co-moving volume the number of uncorrelated regions 
will not be affected by expansion. Now consider the effect of inter-commutings 
(whether or not a loop is produced). Their effect on N can easily be obtained by 
multiplying the inter-commuting rate (obtained from the original analytical model 
Q] ) by the number of regions created by each inter-commuting—this will obviously 
be a positive term. Laguna and Matzner Jll| l have numerically shown that whenever 
two current-carrying strings cross, they inter-commute and a region of intermediate 
current is created. This means that inter-commutings will in general create four 
new regions (see figure [I] (a)). An exception to this is that when regions with size of 
order £ or smaller self-intersect it is possible (see figure [j] (b-c)) that no new regions 
are produced. Also, when the inter-commuting does produce a loop, the regions in 
the corresponding segment are removed from the network, together with one of the 
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newly created ‘intermediate’ regions, so that there will be a negative term in the 
evolution of N (with a similar correction factor at small scales). 

Finally, the only non-trivial issue is that of the dynamics of the currents them¬ 
selves. The simulations of Laguna and Matzner JTj] show that as the result of 
inter-commutings charges pile up at current discontinuities and move with the kinks, 
but their strength decreases with time. Clearly, this indicates that some kind of 
‘equilibration’ process should act between neighbouring current regions, which will 
counteract the creation of new regions by inter-commuting (and help the removal of 
regions by loops). Furthermore, Austin, Copeland and Kibble have shown [|l2j that 
in an expanding universe correlations between left- and right-moving modes develop 
due both to stretching and inter-commuting (particularly when loops form). We 
model this term by assuming that after each Hubble time, a fraction / of the N 
regions existing at its start will have equilibrated with one of its neighbours, 


dN\ 

dt J 


= —fHN ; 


dynamics 


(2.4) 


note that new regions are obviously created by inter-commuting during the Hubble 
time in question, so that / can be larger than unity. Alternatively we can say that 
for a given /, the number of regions that were present in a given volume at a time 
t will have disappeared due to equilibration at a time t + ( fH ) _1 . Presumably the 
only way to find out what / is is by means on numerical simulations (in particular, 
we would expect it to be a model-dependent quantity), although some physical 
arguments can be used to constrain it || . 

However, setting this issue aside for the time being, we obtain the following 
evolution equation for N 


dN =r [P | VooV_ 
dt \£J a L 4 


fHN, 


(2.5) 


where the ‘correction factor’ G has the form 


G = { 2 C (« +2 ) ’ « > 1 

\U { 2(l-2c)a + (2-3c-2a + 4ac)|, f < 1 


( 2 . 6 ) 


loops are assumed to form with a size £(t) = a(t)L{t), with a ~ 1 while the string 
network is is the friction-dominated epoch and a = a sc <C 1 once it has reached 
the linear scaling regime (see [||). Note that when £ £ the net effect of inter- 

commuting and loop production is to remove uncorrelated regions (because each 
loop formed removes a large number of them); otherwise, the net effect is to create 
new regions. 

For what follows it is more convenient to introduce Nl , defined to be the number 
of uncorrelated current regions per long-string correlation length, 



in terms of Nl, (245) has the form 


dN L 

dt 


(3t& - f)HN L + ~-N L + (^G(a7V L ) + | cNj)j V -f-, 


(2.7) 


( 2 . 8 ) 


(with the obvious definition for G{aNif)) where £{ is the relavant friction lengthscale. 
Note that to obtain this equation one needs to use the evolution equation for the 
long-string correlation length L. 
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We should also say at this stage that once the network leaves the friction- 
dominated regime and strings become relativistic other mechanisms (notably ra¬ 
diation) can cause charge losses. Thus we do not expect our toy model to provide 
quantitatively correct answers in this regime, but we do expect it to provide reliable 
order-of-magnitude estimates. 


III. EQUILIBRATION AND SCALING 


We start by d iscussing the case of neutral currents, in which if is given by 
Analysis of (IL8) together with the evolution equations for the long-string network 
reveals two types of behaviour. Firstly, if / is small (that is, equilibration is ineffec¬ 
tive) then Nl grows without limit. In the particular case / = 0, Nl oc t once the 
long-string network has reached the linear scaling regime, meaning that £ oc const. 
On the other hand, for 


/ > fmin = 3 - —— - 1.88 (3.1) 

k + c 

(where we have used the scaling values of c and k obtained from numerical sim¬ 
ulations in the radiation era—see |^] for a discussion) the late time behaviour is 
Nl oc const., £ oc t, that is the superconducting correlation length is scaling lin¬ 
early just like the long-string correlation length; the scaling value of Nl has a fairly 
weak dependence on /, with smaller Nl’s corresponding to larger /’s as expected. 

Note that the fact that correlations cannot obviously be established faster than 
the speed of light (that is, we must have £ < t), means that there is a maximum 
value allowed for /, which at late times (in the radiation era) can be written 

/ < fmax = 3 + kl/ ^ k + g) } 3/2 - 22.4 . (3.2) 


IV. TWO EXAMPLES 

In figure || we plot the result of the numerical integration of (^^) for GUT string 
networks in the friction-dominated epoch, for initial conditions representative of 
string-forming and superconducting phase transitions of first and second order, for 
the cases / = 0 and / = 3, which are representative of the cases where equilibration 
is ineffective and effective, respectively. We are assuming that both of these phase 
transitions occur at around the same (GUT) energy scale. It was also assumed that 
the value of a in the linear scaling regime is a sc ~ 10 -3 (see Martins and Shellard 

§)• 

The differences between the two cases are considerable. If there is no equilibration 
mechanism (/ = 0), £ is conformally stretched during the stretching regime (just like 
the long-string correlation length, L), and so Nl is approximately constant. How¬ 
ever, as inter-commutings start creating new regions Nl begins to increase, growing 
as t 3//2 during the Kibble regime and eventually (once the network reaches the linear 
scaling regime) ends up growing as Nl oc t , which corresponds to £ oc const. As 
expected, in this case the network keeps a ‘memory’ of its initial conditions. On the 
other hand, if there is an equilibration mechanism (/ = 3) then Nl decreases while 
the network is being conformally stretched. In the Kibble regime, the increased 
number of inter-commutings again drives Nl up, and after a has evolved into its 
linear regime value £ itself reaches a scaling value and hence Nl becomes a constant. 
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V. PLASMA EFFECTS 


If the cosmic strings interact with a plasma there will be a further damping term, 
with the corresponding friction lengthscale being given by 


>. = a 

3 p'/ 2 j * 


(5.1) 


where p is the plasma density and J is the string current; at a scale r this is given 
by 


A 1 / 2 

J(r) = e—-— , (5.2) 

r 

N r being the number of uncorrelated current regions at that scale. 

We must therefore determine which of the two damping lengthscales is dominant 
(that is, smaller). It is straightforward to find that at a time t = xt c (t c being the 
epoch of string formation) 


/32tt\ 1/2 

V 3 ) deN^iGp,) 1 / 2 ’ 


where L = jt and d = 47r(7rA/’/45) 1 ^ 2 , (A f being the number of effectively massless 
degrees of freedom). We therefore require this ratio to be less than unity for plasma 
effects to be dominant. 

Now, what initial currents do we expect at t c l If plasma effects are to be impor¬ 
tant, < £ < t at t c : note that £j depends on £, but we can still solve consistently 
for the minimum allowed £j, which will produce the largest current. Using these 
bounds we find in the GUT case 


4 < 


< 60, 


(5.4) 


so that plasma effects are initially subdominat; the lower bound is Gp,-independent, 
whereas the upper one varies as (G/i) -1 / 2 . However, note that the lower bound 
depends on A/", with a ratio of four corresponding to a minimal GUT model A/" ~ 
106.75; one would need A f ~ 10 4 for plasma effects to become dominant—but even 
so, this would only happen in the ‘extreme’ case of a strongly first order string¬ 
forming phase transition and a second order superconducting phase transition. 

Th us w e can in general study the initial stages of superconducting string evolution 
with (T8) and the Aharonov-Bohm frict ion lenghtscale. But then the evolution of 
the currents is such that the ratio fl5.3| ) is at most approximately constant—this 
happens in the ‘extreme’ case of no equilibration, / = 0, in which Nl respectively 
behaves as Nl ex const., Nl ex t 3 / 2 and Nl ex t, in the stretching, Kibble and 
linear scaling regimes. If there is equilibration (that is / ^ 0) then Nl grows more 
slowly than the above and the plasma damping effect loses importance relative to 
Aharonov-Bohm scattering (that is, £j grows faster than £{). 

Hence we find that plasma effects are in general subdominant. However, a way 
of avoiding these constraints is to have further current-building mechanisms, such 
as background magnetic field -this is in fact commonly assumed (and not very well 
explained) in most of the existing litterature (l^] . What we have shown is that if no 
such charging mechanism exists than the currents on the strings are much smaller. 
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VI. CONCLUSIONS 


In this paper we have constructed a simple ‘toy model’ for the evolution of currents 
on cosmic strings, based on the dynamics of a ‘superconducting correlation length. 
The different dynamical processes affecting it were considered, and a simple study 
of the solutions of the model has revealed that the process of equilibration is crucial 
in the evolution of the currents. We have found that a fairly efficient equilibration 
mechanism is needed for the superconducting correlation length to scale. Whether 
or not such a mechanism is available in realistic model is an issue which cannot be 
addressed at the moment—a detailed numerical investigation will clearly be needed 
for that. 

We have also studied the importance of plasma damping effects on the strings. 
We find that these are typically negligible when compared to the ‘usual’ Aharonov- 
Bohm scattering. However, an exception to this might happen in the presence of 
background magnetic fields, either of ‘primordial’ origin or generated (by a dynamo 
mechanism, say) after protogalaxies start to collapse (an example of the later is 
the Chudnivsky-Vilenkin scenario [^3|). Even in this case, for typical values of 
cosmological or astrophysical magnetic fields, the two damping lengthscales are not 
too different. 

While this paper was being prepared, a rather less accurate discussion of plasma 
effects on string dynamics appeared in the litterature p|. This is in broad agree¬ 
ment with our results if one allows for the fact that in [ no consideration is given 
to the dynamics of currents—in fact current is simply treated as a constant. 

Our extended quantitative model allows reliable estimates of the currents on 
strings to be made at all times. Thus a more detailed analysis of some of the 
cosmological scenarios involving superconducting cosmic strings becomes possible. 
An outstanding example the question of the abundance of vorton remnants—we 
will address it in a forthoming publication Q. 
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FIG. 1. Some relevant inter-commuting configurations. The arrows mark the limits of 
regions with correlated currents. Plot (a) shows a typical inter-commuting creating four 
new current regions, while (b-c) show than on scales smaller than the current correlation 
length loop production may (c) or may not (b) remove current regions from the long-string 
network. 
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FIG. 2. The evolution of the number of uncorrelated current regions per long-string 
correlation length, Nl, for the cases / = 0 (top) and / = 3 (bottom) assuming that the 
orders of the string-forming and superconducting phase transitions are respectively: 1st 
& 1st (solid lines), 1st & 2nd (long dashed), 2nd & 1st (short dashed) and 2nd & 2nd 


(dotted). Time is in orders of magnitude from the epoch of string formation. 
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